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We study theoretically the effects of anisotropy on the thermoelectric performance of p-type AgBiSe 2 . We 
present an apparent realization of the thermoelectric benefits of one-dimensional “plate-like” carrier pocket 
anisotropy in the valence band of this material. Based on first principles calculations we find a substantial 
anisotropy in the electronic structure, likely favorable for thermoelectric performance, in the valence bands of 
the hexagonal phase of the silver chalcogenide thermoelectric AgBiSe 2 , while the conduction bands are more 
isotropic, and in our experiments do not attain high performance. AgBiSe 2 has already exhibited a ZT value of 
1.5 in a high-temperature disordered fee phase, but room-temperature performance has not been demonstrated. 
We develop a theory for the ability of anisotropy to decouple the density-of-states and conductivity effective 
masses, pointing out the influence of this effect in the high performance thermoelectrics Bi 2 Te 3 and PbTe. From 
our first principles and Boltzmann transport calculations we estimate the performance of p-type AgBiSe 2 . 


I. INTRODUCTION 

Anisotropy is a substantial contributor to many phenomena 
of technological importance in condensed matter physics. For 
example, today’s high performance magnets owe their perfor¬ 
mance in large part to a sizable magnetocrystalline anisotropy, 
deriving from spin-orbit coupling, the associated crystal field, 
and a non-cubic crystal symmetry. 

Anisotropy is also important for thermoelectric perfor¬ 
mance, or the figure-of-merit ZT (defined in Appendix A). 
Very recently! Mecholsky et al found that band “warping”, or 
a non-analyticity of the effective mass tensor near the band 
edge, can have a substantial effect on the electronic transport 
relevant for thermoelectric behavior. In that work as well as 
the recent review by Shakouri! it was pointed out that the 
existence of a band edge with different masses in different 
directions favors thermoelectric performance. This was also 
shown in Ref. [H and the benefits of anisotropy to thermo¬ 
electric transport, with r^ards to FeAs 2 and LiRh 204 , were 
shown in Refs. 0 and 0 In fact this basic scenario, ap¬ 
plied to artificial structures, had been foreseen by Hicks and 
Dresselhaus^ in works on thermoelectric transport in super¬ 
lattices and nanowires. It was also shown^ that one can realize 
the beneficial effects of low-dimensional electronic structures 
even in isotropic materials. 

Numerous authors have noted^rJ? the effect of carrier 
pocket degeneracy on thermoelectric performance, noting that 
such degeneracy improves electrical conductivity without sac¬ 
rificing the Seebeck coefficient, both indispensable compo¬ 
nents of a useful thermoelectric. Such degeneracy is nearly 
universally rooted in the placement of a band extremum away 
from the F point since the crystal symmetry then dictates the 
appearance of multiple carrier pockets. This crystal symme¬ 
try also ensures that anisotropy in the electronic transport re¬ 
sulting from a single carrier pocket, which we term “pocket 
anisotropy”, does not sacrifice the overall electronic transport 
of the crystal. Examples of such pocket anisotropy include the 
valence bands of the high performance thermoelectrics Bi 2 Te 3 
and PbTe. The valence bands of both these materials contain 
band edges away from the F point with threo^S (for Bi 2 Te 3 ) 
and two!! (for PbTe) disparate effective masses in different di¬ 
rections. However, the transport is isotropic in the basal plane 


for Bi 2 Te 3 and for all three dimensions in cubic PbTe. Our 
calculations described below find essentially isotropic over¬ 
all electrical transport for p-type AgBiSe 2 despite substantial 
carrier pocket anisotropy. 

In addition, the benefits of mass anisotropy increase as the 
dimensionality of the electronic structure decreases. For ex¬ 
ample, an isotropic three dimensional parabolic band has a 
v/E dependence of the density-of-states N{E) on energy E, a 
two-dimensional parabolic band has a step-function behavior, 
and a fully one-dimensional band has an behavior, im¬ 

plying a diverging density-of-states (DOS) at the band edge. 
Note that a one-dimensional band takes on a “sheet” or “plate¬ 
like” isoenergy or Fermi surface. While in real materials such 
a DOS divergence does not occur due to a lack of complete 
one-dimensionality, it is clear there are particular benefits of 
such a one-dimensional structure. 

In this work we present calculations of such one¬ 
dimensional “plate-like” carrier pocket anisotropy in the 
valence band of the silver chalcogenide semiconductor 
AgBiSe 2 . This material crystallizes in a hexagonal structure 
at room temperature, with successive transitions to a rhom- 
bohedral structure at 420 K and to a disordered fee structure 
at about 600 K. Recent work!!"— shows a ZT value for n- 
type AgBiSe 2 of 1.5 just above the cubic-rhombohedral phase 
transition. We find from theory that good room-temperature 
p-type thermoelectric performance may be possible, estimat¬ 
ing a p-type ZT value of 0.4-0.7, considering only electronic 
optimization. 

Ordinarily one-dimensional electronic structures are asso¬ 
ciated with electronic instabilities such as charge or spin- 
density waves. These instabilities typically lead to insulat¬ 
ing ground state behavior. Here, by contrast, we have a one¬ 
dimensional electronic structur&i associated with a metallic 
ground state, that we will show to be highly beneficial for 
thermoelectric performance. 

We observe that very few materials have ZT ^ 1 in the 
important temperature range between 200 and 350 K, where 
many heating and cooling applications exist. Substantial ef¬ 
forts to find a substitute for the prototypical thermoelectric 
Bi 2 Te 3 , which shows a ZT of unity around room tempera¬ 
ture, have been made for nearly 50 years. Hence the find¬ 
ing of a candidate for such performance levels, assuming full 
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optimization, in this temperature range is of substantial im¬ 
portance. This remains true even despite the difficulty, de¬ 
scribed below, of doping into the valence band of AgBiSe 2 . 
We note that both Bi 2 Te 3 and AgBiSe 2 contain expensive el¬ 
ements (Tellurium and Silver). However, Silver is more avail¬ 
able, with worldwide annual production is 26,000 tons, com¬ 
pared to Tellurium, with annual production of 100 tons^. This 
parallels the 75 times larger mass abundance^ in the Earth’s 
crust of Silver relative to Tellurium. 

The remainder of this paper is organized as follows: in Sec¬ 
tion II, we present our experimental procedure; in Section III 
we present our experimental data, along with a comparison of 
the calculated n-type and p-type electronic structures; in Sec¬ 
tion IV, we present a more general consideration of anisotropy, 
followed in Section V with our first principles theoretical re¬ 
sults and our estimate of ZT, and our conclusions in Section 
VI. Details on the estimation of ZT are presented in Appendix 

A, and a demonstration of the isotropy of the Seebeck coeffi¬ 
cient within an effective mass approach is given in Appendix 

B. 


II. EXPERIMENTAL PROCEDURE AND EFFORTS TO 
ATTAIN p-TYPE DOPING 

AgBiSe 2 samples were prepared by reacting high-purity el¬ 
ements in evacuated silica ampoules. For the ‘as-grown’ sam¬ 
ple, solidihcation occurred slowly (1 deg/hr) to promote large 
grains. The sample was then cooled at 5deg/hr between 310°C 
and 260°C, followed by a 48h anneal at 260°C to minimize 
defects and reduce residual stresses. Due to the formation of 
voids during growth and and cracks during processing, only 
relatively small pieces were obtained. These difficulties lead 
to the synthesis of polycrystalline samples, which also allow 
for a more rapid screening of potential dopants. Following re¬ 
action in the melt, samples were ground in a He glove box and 
hot-pressed in a graphite die at temperatures near 450°C with 
a pressure of approximately 10,000 psi. This process resulted 
in samples with geometric densities of approximately 95% of 
the theoretical density. 

The Seebeck coefficient was used as a screening tool to ex¬ 
amine the influence of the various substitutions and process¬ 
ing conditions. The Seebeck and Hall coefficients were mea¬ 
sured in a Quantum Design Physical Property Measurement 
System, using the Thermal Transport and Resistivity Options, 
respectively. Thermal measurements were performed using 
gold-coated copper leads attached to the sample with H20E 
Epo-Tek silver epoxy. Hall data were obtained using a stan¬ 
dard four-wire conhguration, with 0.0508 mm Pt wires spot 
welded to the sample and maximum fields of ± 6 T were em¬ 
ployed; electrical resistivity was collected during the same 
measurement as the Seebeck coefficient and thermal conduc¬ 
tivity. Scanning electron microscopy (SEM) and energy dis¬ 
persive spectroscopy (EDS) measurements were performed in 
a Hitachi TM-3000 microscope equipped with a Bruker Quan- 
tax 70 EDS system. 

In order to manipulate the electrical properties, polycrys¬ 
talline samples containing S, Te, and Pb were made. Sim¬ 


ilarly, a sample of nominal composition Agi iBio. 9 Se 2 was 
produced from the melt to examine the potential for instrin- 
sic hole doping via bismuth deficiency. SEM coupled with 
EDS revealed a Ag-rich phase at apparent grain boundaries in 
this Agi iBio. 9 Se 2 sample, which indicates the tendency of the 
phase to form near the stoichiometric composition AgBiSe 2 . 
This sample was not considered further due to our desire to 
probe, as much as possible, only the intrinsic properties of 
AgBiSe 2 . 

In our experimental results, AgBiSe 2 naturally forms n- 
type. As-grown AgBiSe 2 was found to have a Hall carrier 
density iih near 2.3 x lO'^cm^^ at room temperature, and little 
temperature variation was observed down to 25K. At 300 K, 
this as-grown sample had a Seebeck coefficient of -138 /jV/K 
and a Hall mobility of /j//=45cm^A^/s. Our polycrystalline 
AgBiSe 2 had a Seebeck coefficient of Ri-400/jV/K at 300 K, 
and an electrical resistivity that increased with decreasing 
temperature. These results suggest that our polycrystalline 
AgBiSe 2 is very lightly doped, or nearly an intrinsic semi¬ 
conductor. Thus, AgBiSe 2 can be formed near the insulating 
limit, which makes obtaining p-type conduction more plausi¬ 
ble. 

Pb has previously been showr^ to induce p-type conduc¬ 
tion in AgBiSe 2 . However, our Pb-doped sample of nominal 
composition AgBio, 95 Pbo.o 5 Se 2 had a Seebeck coefficient of 
~-95pV/K, suggesting a higher electron concentration than 
the as-grown sample or the polycrystalline AgBiSe 2 . The 
apparent increase in free electrons with a p-type dopant is 
surprising, and may be related to inhomogeneities within the 
sample. However, the absence of p-type conduction at low Pb 
concentrations was reported by Pan et al2^, though they ob¬ 
served a cross-over to p-type conduction near x = 0.02. We 
did not continue to investigate Pb-doping after this initial hnd- 
ing. Instead, we briefly examined the influence of S and Te 
substitutions. 

The n-type behavior of AgBiSe 2 may be caused by Se 
vacancies. To gain insight into this potential mechanism, 
we considered S and Te substitutions for Se. Interestingly, 
our sample of nominal composition AgBiSei 8 Teo .2 had a 
room temperature Seebeck coefficient of -475pV/K and re- 
sitivity that increased with decreasing T. Compared to the 
-400pV/K observed for the undoped, polycrystalline sample, 
this would suggest that Te substitution drives the systems to¬ 
wards a charge balanced state due to its lower vapor pres¬ 
sure (fewer anion vacancies). However, we remain cautious 
because AgBiTe 2 is also naturally n-typd22.. Also, Seebeck 
coefficient measurements become more difficult and absolute 
errors increase for resistive samples. Yet, sulfur substitution 
pushed the system in the opposite direction (more metallic), 
consistent with a relative increase in elemental vapor pressure. 
Unfortunately, there was no trend in these data and the more 
heavily substituted sample had a larger 300 K n-type Seebeck 
coefficient (-87pV/K for 5% S, and -175pV/K for 10%S). As 
such, we cannot draw any conclusions regarding the sulfur 
substitution, though we remain optimistic about coupling Te 
substitution with non-isoelectronic dopants. 

A detailed study of the defect chemistry in these materi¬ 
als could provide the insights required to achieve the desired 


3 


levels of /?-type conduction. A more complete picture would 
require additional experimental and theoretical efforts, con¬ 
sidering the variety of compositions in this family of I-V-VI 2 
chalcogenides. 

The large Seebeck coefficients and semiconducting behav¬ 
ior of the resistivity of our polycrystalline samples thus sug¬ 
gest that an adequate p-Wpe dopant might be found. We also 
note that in Refs. |23|and0/7-type behavior - i.e. positive See¬ 
beck coefficient - was observed in the low-temperature hexag¬ 
onal phase, so our results should not be taken to imply that 
p-type doping is impossible to achieve in AgBiSe 2 , merely 
that more involved effort will likely be required to achieve 
this. In the literature there are numerous semiconductors such 
as CrSi 2 and Bi 2 Se 3 ^^'^^ which were originally found to ex¬ 
hibit a strong doping type preference, which were later found 
possible to dope both n-type and p-type under the proper cir¬ 
cumstances. Since AgBiSe 2 has not undergone much exper¬ 
imental study to date, we think it likely that the difficulty in 
attaining p-type doping will be successfully addressed by fu¬ 
ture efforts. We note finally that the relatively narrow calcu¬ 
lated band gap of 0.52 eV (see Section IV) argues in favor of 
the likelihood of attaining both doping types. 


III. EXPERIMENTAL DATA AND THEORETICAL 
COMPARISON OF «-TYPE AND p-TYPE 

In Figure 1 we show the T-dependent Seebeck coefficient, 
resistivity and thermal conductivity for the as-grown sam¬ 
ple, along with (left panel) our hrst principles calculation 
of the Seebeck coefficient for this sample. The measured 
Hall number, -0.27 cm^/Coul, would correspond to an elec¬ 
tron concentration n of 2.3 x lO'® cm^^ in the case of an 
isotropic parabolic band, although Fig. 2 (top) depicts a con¬ 
duction band electronic structure signihcantly differing from 
isotropic, so that the chemical and Hall-inferred carrier con¬ 
centrations may differ. We hnd a good ht at a chemical carrier 
concentration n of 2.52 x 10*® cm^^ indicating the accuracy 
of our theoretical approach, which uses the constant scattering 
time approximation^. 

The right hand panel of Fig. 1 shows the thermal conduc¬ 
tivity and resistivity measurements. We see that the thermal 
conductivity K of AgBiSe 2 is extremely low - it reaches a min¬ 
imum value of 1.1 W/m-K around 200 K, indicating low lat¬ 
tice thermal conductivity - a key parameter of a useful thermo¬ 
electric. In fact, using the resistivity data and the Wiedemann- 
Franz relationship we hnd the lattice term at 200 K to be just 
1.0 W/m-K, which is lower than that of the better-known high 
performance thermoelectrics such as Bi 2 Te 3 and PbTe. This 
value is slightly higher than that found by Nielsen et al for 
rhombohedral AgBiSe 2 r^, presumably due to the difference in 
physical structure. The slight upturn of K above 200 K likely 
rehects radiative effects not corrected for here. 

The resistivity data shows the Fermi-liquid T^ dependence 
at the very lowest temperatures below 50 K, crossing over to 
T-linear behavior for temperatures above 100 K. This T-linear 
behavior is characteristic of electron-phonon scattering and al¬ 
lows us to estimate the electron phonon coupling constant A. 



FIG. 1. (Left panel) The measured and calculated thermopower of 
AgBiSe 2 ; right panel, the measured resistivity (blue triangles) and 
thermal conductivity (red squares). 


for this sample. We will later use this to estimate the perfor¬ 
mance of p-type AgBiSe 2 . From our theoretical calculations 
we find the plasma frequency squared ( 0 ^ at the modeled dop¬ 
ing level to be 0.096 eV^. When we combine this with our re¬ 
sistivity data and the theoretical relationship connecting A, ®/>, 
and the resistivity from Ref. we hnd an electron-phonon 
coupling constant of 0.49, and an associated electronic scat¬ 
tering time of 8 x 10^*^ s. 

Although we were able to able to obtain signihcant n-type 
doping levels, our data suggest that in this temperature range 
n-type AgBiSe 2 is not likely to be a high performance thermo¬ 
electric. The 300 K power factor for our as-grown sam¬ 
ple is just 0.3 mW/m-K^, or less than 10 percent of the value 
of optimized Bi 2 Te 2 . Although one might achieve some gain 
in the power factor in a more lightly doped sample, this will 
not likely raise 300 K ZT substantially from the ~ 0.1 value 
achieved here (a rough estimate hnds optimized ZT values of 
less than 0.2). 

It is of interest to understand the reason for this. Pre¬ 
sented in Figure 2 (top) is a plot of the first-principles cal¬ 
culated isoenergy surface for n-type AgBiSe 2 for a doping 
n = 1.23 X 10^® cm^^. The plot depicts a cylindrical body 
whose width is quite comparable to its height. While the 
cylindrical surface is suggestive of two-dimensionality and 
hence mass anisotropy, the shape does not lend itself to a large 
surface to volume ratio, which we have argued elsewhere^ to 
be favorable for high ZT. The structure is also substantially 
lacking in degeneracy (technically there is a two-fold degener¬ 
acy as the band edge is at the A point). Given these factors and 
the lack of the favorable complexity described in Ref. [sH it is 
perhaps not surprising that the performance levels of hexago¬ 
nal n-type AgBiSe 2 are comparatively low. 

The situation is rather different for p-type, however. Figure 
2 (bottom) presents a plot of the isoenergy surface for p-type 
AgBiSe 2 at the same carrier concentration as depicted for n- 
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FIG. 2. The calculated Fermi surface of hexagonal AgBiSe 2 at 
an electron doping n = 1.23x 10^®cm^^) (top) and the same carrier 
concentration for p-type (bottom). The bottom plot is rotated slightly 
to better depict the degeneracy and mass anisotropy. 


type. This doping is likely near the optimal p-type doping 
for room-temperature ZT. The Fermi surface consists of six 
platelike structures that exhibit the feature mentioned in the 
Introduction: a perpendicular mass much less than the parallel 
mass. In fact, in the plane of the “plate”, the effective mass is 
approximately 5 times the perpendicular mass. There are ad¬ 
ditionally highly elongated features which stretch out to a sec¬ 
ondary band maximum at a non-symmetry point on the zone 
surface. We also note the six-fold degeneracy, a combination 
of the hexagonal symmetry and the position of the VBM away 
from r, and away from high-symmetry zone boundary posi¬ 
tions. Finally, one observes the orientation of the “plates”: 
they are located at an angle to the c-axis, not parallel to it, so 
that the favorable transport properties of this feature will ex¬ 
tend to the c-axis transport as well. We will see that p-type 
transport in this material is expected to be virtually isotropic. 

In the next section we give a general theoretical descrip¬ 
tion of the benefits of carrier pocket anisotropy in enhancing 
thermoelectric performance, as appear to be active in p-type 
AgBiSe 2 , specifically the role of such pocket anisotropy in de¬ 
coupling the density-of-states or thermopower effective mass 
from the conductivity effective mass. 


IV. THEORY OF ROLE OF POCKET ANISOTROPY IN 
DECOUPLING DENSITY-OF-STATES AND 
CONDUCTIVITY EFFECTIVE MASSES 


Here we show that any semiconductor band edge pocket 
anisotropy that can be expressed in terms of an ellipsoidal 
effective mass tensor has substantial benefits to thermoelec¬ 
tric performance, provided that the band edge is not located 
at the F point and the crystal obeys a certain minimum sym¬ 
metry. Larger anisotropies, in addition to the effects of band 
degeneracy^., are shown to be more beneficial in this regard. 

We begin with the canonical expressions for the temper¬ 


ature and chemical-potential-dependent thermopower tensor 
‘^aa(T,/j) and the electrical conduct ivit y tensor Oaa 
which we reproduce here from Ref. [3^ For simplicity we 
assume that both are diagonal in the spatial indices a and also 
assume only one band contributes to transport. Then we have: 


Saa{T,/j) 


Vaa(T,/j) 


( 1 ) 


with 


OaaiT,^) = -IdE ( 2 ) 

and 

Va.a(T,A<) = ^/ Oaa{E){E - ^i)^-^^^^dE (3) 

Here am p(£’) is the transport function which is written as 

Oaa{E) = „ ,T,'k§(£' - (4) 


and / is the Fermi function. Note that for an anisotropic 
parabolic band, in which the relaxation time x depends on k 
through F; k, the above expression (Eq. 4) for the transport 
function can be written as 


Oa,a(£) = e\l{E)%a{E)N{E) (5) 


Here the effective mass anisotropy is incorporated into the di¬ 
rectional index a. 

As described in Appendix B, this effective mass anisotropy 
nonetheless yields an isotropic Seebeck coefficient, so long as 
the bands are taken as parabolic, only carriers of one sign con¬ 
tribute to transport, and the scattering time is taken as depend¬ 
ing only on energy, not on direction. We note also that, though 
difficult to obtain in practice, a highly anisotropic Seebeck co¬ 
efficient is of interest in the study of transverse thermoelectric 
effects^. 

We now analyze this problem through the concept of effec¬ 
tive mass. From standard references^ the electrical conduc¬ 
tivity tensor may be written as 


o = 


E d k 1 

47t3 


( 6 ) 


where is an effective mass tensor and the sum is taken 
over occupied levels. The standard electronic conductivity 
formula (where we have retained the directional indices) can 
then be written as 


eta,a — 


ne^Xg 

mg.a 


(7) 


and contains a conductivity effective mass nig g to which the 
conductivity is inversely proportional. For high electrical con¬ 
ductivity (a prerequisite for good thermoelectric performance) 
one therefore desires small effective masses, irrespective of 
the scattering time x,. 

The situation is very different for the thermopower S. In the 
fully degenerate limit the Mott formula gives the thermopower 
as 


Sa,a{n,T) — —- hgT 

3 e 


ks (d\og{(3aa(E)) 


dE 


\E=Ef 


( 8 ) 
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For an anisotropic parabolic band, one has (from Appendix 
B) that 


<yaa{E) = (9) 


where E is measured from the band edge and g 3 is a function 
of the band masses. The logarithmic derivative in the above 
equation is simply 3/2Ef. Ep may be rewritten in terms of 
the carrier concentration n and effective masses nix ,and ni^ 
by noting that n can be written as (using the formula for the 
volume of an ellipsoid and including a factor of 2 for spin 
degeneracy) 


n = 


2 47t 
(27t)3 T 


^F,x^F,y^F,z 


( 10 ) 


where hkp.i = y/lniiEp with m, the effective mass in direction 
/, yielding 


n = 


y/2,mxmym^Ey^ 


( 11 ) 


Combining the above equations and rewriting {nixniym^^/^ as 
niDos, we finally have 


S{T,n) 


eh^ 


moosT 



( 12 ) 


Eor the thermopower larger effective DOS masses are benefi¬ 
cial, which runs completely counter to the benefit to the con¬ 
ductivity of smaller masses. This is one major fundamental 
problem that must be overcome to achieve high thermoelec¬ 
tric performance. 

However, as suggested by the notation, there is no particu¬ 
lar reason that nia and moos must be equal. They are in fact 
equal only for a single isotropic parabolic band. This distinc¬ 
tion is well-known for the case of band degeneracy, where a 
band edge of degeneracy N enhances^ moos by a factor 
without affecting Here we explore in addition the effect 
of anisotropy on the relation of these two masses. We find that 
anisotropy allows for small conductivity effective masses iria, 
heightening the electronic conductivity, but large density-of- 
states effective masses moos, enhancing the thermopower. 

One can see the reason for this quite easily. Eor simplic¬ 
ity, (the effect is similar for three distinct masses) consider 
an anisotropic parabolic band edge represented by a three- 
dimensional ellipsoid of revolution, characterized by two ef¬ 
fective masses, a radial mass m[| and a longitudinal mass m±. 
The DOS effective mass for this situation is given by moos 
= (m|mx)'/^. The conductivity effective masses in each of 
the two directions are given by m\\ and m±. However, for 
a crystal in which the band edge is not at the F, or similarly 
low-degeneracy, point, there will generally be some band edge 
degeneracy which respects the crystal symmetry. As depicted 
in Ref. this means that even an anisotropic band edge can 
result in isotropic transport in a cubic material. The same gen¬ 
eral idea applies for planar transport for a layered material 
such as Bi 2 Te 3 . One then sees that the (isotropic) conductiv¬ 
ity effective mass is given by 


3/mo = 2/m||-b l/m^ (13) 


»*DOS 


/Wo- 



FIG. 3. The ratio nioos/nia- Here r = . 


so that 


3m[|mx 
2mx -\-m 


(14) 


As is clear from the expression, this conductivity mass m^ 
has a rather different dependence on nip and m[| than the 
density-of-states mass m^^^. Indeed, the ratio moos/fn^, 
which is effectively a figure-of-merit for the effects of elec¬ 
tronic anisotropy, is simply 


moosl'^G 


2mp -b m|| 

3(m[|m^) 

(2/3)r^/^-b(l/3)r^2/^ 


(15) 

(16) 


In the last expression we have expressed mj^Q^/m^ in terms 
of r where r is the mass ratio mx/m||. Now, this nioos/'^a 
ratio can be substantially different from unity if m|| and m± 
are very different. For a case where nip — 24m[| (a value 
present in the conduction band of GeTe^, this ratio approaches 
the value 2. Additionally, in the opposite limit, in which the 
band edge forms a flattened, highly prolate, rather than oblate, 
spheroid, still larger effects are present. Figure 3 depicts a plot 
of the DOS to conductivity mass ratio as a function of the ratio 
r — mx/m[[. For this ratio = 1/100 one finds a DOS to con¬ 
ductivity mass ratio of nearly 10, which is clearly beneficial 
in enhancing both thermopower and electrical conductivity. 

As examples of these concepts, we consider the p-type band 
masses in the high performance thermoelectrics Bi 2 Te 3 and 
PbTe. The band edge in p-type Bi 2 Te 3 has three distinct 
masses^S, which take the values of 0.064 mo, 0.196 mo, and 
0.73 mo, where mo is the free electron mass. These yield a 
band edge moos of 0.209 mo and m^ of 0.136^ so that their 
ratio is over 50 percent enhanced relative to an isotropic band 
edge. When one includes the six-fold degeneracy of the band 
edge in the DOS mass the effective mass ratio becomes 3.87, 
which is surely a major contributor to the performance of this 
material. Similarly, the L-point band edge in PbTe contains 
a radial mass^ of 0.022 mo and longitudinal mass 0.31 mo. 
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yielding an nioos/tnc, of 1.66, or 4.19 if the fourfold L-point 
degeneracy is included, which likely contributes to the excep¬ 
tional thermoelectric performance^ of this material. 

V. POSSIBLE REALIZATION OE BENEEICIAL EEEECTS 
OE CARRIER POCKET ANISOTROPY: p-TYPE AGBISE 2 

We focus here on the low temperature hexagonal phase of 
AgBiSe 2 , calculating its properties with the first principles 
density functional theory code WIEN2K^, as well as Boltz¬ 
mann transport properties via the Boltztrap code^, within the 
generalized gradient approximation (GGA) of Perdew, Burke 
and Ernzerhofii. Eor these calculations we use a modification 
of the GGA known as a modified Becke Johnson potential^ 
which gives accurate band gaps^^f^ , a matter of great impor¬ 
tance for the transport properties. The linearized augmented 
plane wave (LAPW) basis was used, with LAPW sphere radii 
of 2.5 Bohr for all atoms, and an RK,nax of 9.0 was used. Here 
RKmax is the product of the sphere radius and the largest basis 
wave vector. Approximately 1000 k-points in the full Bril- 
louin zone were used for the self-consistent calculations run 
to convergence and approximately 10,000 points were used 
for the transport calculations. The internal coordinates were 
relaxed using the standard GGA until forces were less than 2 
mRyd/bohr. All other calculations, the relaxations excepted, 
included spin-orbit. 



FIG. 4. The calculated band structure of hexagonal AgBiSe 2 . 

Eigure 4 depicts the calculated band structure of hexagonal 
AgBiSe 2 . Erom the figure, we see that this material is a semi¬ 
conductor of band gap 0.52 eV, with conduction band mini¬ 
mum at the A point and a valence band maximum located be¬ 
tween the r and K points (the Brillouin zone for these points is 
found in Ref. 0). The valence band, in addition contains sub¬ 
stantial degeneracy, with at least four subsidiary maxima less 
than 100 meV from the band edge. These maxima all origi¬ 
nate from the same general Eermi surface structure, and in ad¬ 
dition we note a wide range of masses for these maxima, with 
the maximum between T and H having the heaviest mass and 
the adjacent maximum between H and L having the lightest. 


This accords with the discussion in the Introduction regarding 
the anisotropy in band extrema. To demonstrate the potential 
of this material for high ZT, in Figure 5 we present a plot of 
the 300 K calculated power factor vs. carrier concentration 
c/t and compare with p-type Bi 2 Te 3 , the highest performing 
room temperature thermoelectric. Although the p-type ther¬ 
mopower and conductivity of AgBiSe 2 are nearly isotropic, 
Bi 2 Te 3 has substantial anisotropy, as is well known. To en¬ 
sure a fair comparison we present the conductivity-averaged 
quantities (see Ref. IH for details on these quantities) since in 
a typical experiment a polycrystalline sample is used, which 
tends to average the transport over the principal axes. 

The figure depicts 5 ^c/t results much larger than those of 
Bi 2 Te 3 . In particular, at the doping where the thermopower is 
200p V/K (indicated by the asterisks), the of AgBiSe 2 

is double that of Bi 2 Te 3 . This doping, which is p= 2 x lO'® 
cm^^ forBi 2 Te 3 and 1.2 x for AgBiSe 2 , is approxi¬ 

mately the doping of optimized ZT for Bi 2 Te 3 and is probably 
at or near optimal for AgBiSe 2 . While one cannot exclude 



FIG. 5. The calculated power factor divided by scattering time, 
S^c/t of p-type Bi 2 Te 3 and AgBiSe 2 . The asterisks indicate the 
doping where the calculated thermopower is 200 p V/K, which is 
usually near the optimal doping for a high performance thermoelec¬ 
tric. 

scattering time differences from the comparison, this result 
suggests that AgBiSe 2 may show comparable power factors 
to those of Bi 2 Te 3 , and when combined with the observed low 
lattice thermal conductivity of AgBiSe 2 ^ suggests a high po¬ 
tential for room-temperature thermoelectric performance. 

In Figure 6 we depict the calculated thermopower for 
AgBiSe 2 as a function of carrier concentration. Note that the 
p-type thermopower is virtually isotropic. We have indicated 
the range of thermopowers between 200 and 300 pV/K as this 
is the general range of thermopower over which ZT is maxi¬ 
mized. From this we find optimal 300 K p-type doping levels 
of 3 X 10*® — 1.3 X 10^**cm^^. For n-type the optimal doping 
levels are much lower, indicating a lower likelihood of good 
thermoelectric performance. 

In Figure 7 we present the calculated 300 K electrical con- 
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n,p (cm'^) 


FIG. 6. The calculated thermopower of hexagonal AgBiSe 2 at 300 
K. 



FIG. 7. The calculated electrical conductivity ratio Oc/Oab and 
thermopower ratio Sc/Sah of hexagonal AgBiSe 2 at 300 K. 


ductivity and thermopower ratios Oc/Oab and Sc/Sab. In view 
of the great carrier pocket anisotropy depicted in Figure 2 
(bottom), it is of interest that the thermopower and electri¬ 
cal conductivity of p-type AgBiSe 2 vary by no more than 10 
percent from the c-axis to the plane. Note that the electrical 
conductivity, with respect to an unknown scattering time x, 
is calculated as the transport function N{E)vj{E) integrated 
with the derivative of the Fermi function, a function sharply 
peaked around Ep. Here N{E) is the density of states and v, 
the Fermi velocity in direction i. Hence the anisotropy of the 


conductivity is a good measure of the anisotropy of the trans¬ 
port function itself. 

We note that, as suggested by Figures 6 and 7, in general 
the overall global electronic transport anisotropy need not be 
as large as the local carrier pocket anisotropy. Only in the case 
of a single band extremum are these two anisotropies equal. 
When multiple extrema, as depicted in the bottom panel of 
Figure 2, are present the overall anisotropy can be much less 
than the carrier pocket anisotropy, due to the varying relative 
orientations of the band extrema dictated by the crystal sym¬ 
metry. Thus despite the great variation in directional effec¬ 
tive mass in the hole pockets in the bottom of Figure 2, the 
overall electronic transport in p-type AgBiSe 2 is very nearly 
isotropic. 

Actual values of the electrical conductivity depend on the 
electronic scattering time x, which can vary substantially both 
from one material to another and by doping level within a par¬ 
ticular material. As is well known, for temperatures above 
the Debye temperature^ x is generally inversely proportional 
to the dimensionless electron-phonon coupling constant A,, as¬ 
suming the absence of extrinsic factors such as grain boundary 
scattering. While A is not directly available from the first prin¬ 
ciples calculations we have performed, we estimated above its 
value as 0.49 from experimental measurements described in 
Section II. This section also presented measurements of the 
thermal conductivity, which allows estimation of the potential 
ZT of p-type AgBiSe 2 . 

As described in Appendix A, from these data we estimate 
the 300 K ZT value of optimally doped p-type AgBiSe 2 as 
0.4 - 0.7. These values consider only optimization of the 
electronic transport; optimization of the lattice transport can 
be expected to yield additional performance benehts. These 
promising performance values are representative of the posi¬ 
tive effects of carrier pocket electronic anisotropy in produc¬ 
ing favorable thermoelectric behavior. 


VI. CONCLUSION 

Carrier pocket electronic anisotropy is seen to positively 
impact thermoelectric performance, provided that the crystal 
obeys a minimum symmetry and the relevant band edge 
is located away from the F point so that the material can 
experience the benehcial effects of both band degeneracy 
and a large effective mass ratio in the ellipsoidal effective 
mass tensor, p-type AgBiSe 2 appears to be a material in 
which these benefits are present, with “plate-shaped” Fermi 
surfaces, a large predicted power factor and ultimately a 
predicted 300 K ZT value of 0.4 - 0.7, considering only 
optimization of the electronic properties. The favorability 
of this material suggests that other materials with such 
anisotropic features may exist; indeed, two of the best 
performing thermoelectrics Bi 2 Te 3 and PbTe, exhibit a great 
deal of such anisotropy. Searches for other potentially high 
performance thermoelectrics with these anisotropic behaviors 
may therefore be of interest. 





















Appendix A ; Method for estimation of thermoelec¬ 
tric performance of p-type AgBiSe2 
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Thermoelectric performance is measured as 


zr = 


S^oT 


(17) 


K 

Here S is the Seebeck coefficient, K the thermal conductiv¬ 
ity, and c the electrical conductivity. To make quantitative 
estimation of thermoelectric performance we need the values 
of these three quantities. We do this using a combination of 
our first principles calculations, the experimental data just de¬ 
scribed, and a few assumptions. 

Although in p-type AgBiSe 2 the Seebeck coefficient 
and electrical conductivity are calculated to be essentially 
isotropic, we describe here briefly the application of Eq. (17) 
in the anisotropic case. For a single crystal sample, the above 
expression for ZT directly applies for the Seebeck coefficient, 
electrical and thermal conductivity measured in a particular 
direction, leading to a value of ZT which is directionally de¬ 
pendent. The more typical situation, as in our experimental 
work, is that of a polycrystalline sample with random grain 
orientation. In that case ZT of the sample is isotropic, and 
the values of S and a which enter the ZT expression are as 
follows (we assume diagonal S and c tensors); 


^ _ Sxx^xx + ^yy^yy T 

^avg 


^avg — — 


Jyy- 


(18) 

(19) 


Note that the above expression neglects extrinsic effects such 
as grain boundary scattering that may reduce Oavg from the av¬ 
erage of the corresponding single-crystal direction-dependent 
conductivities. A detailed consideration of the effect of trans¬ 
port anisotropy on ZT may be found in Ref. 

Returning to calculating the ZT of p-type AgBiSe 2 , our 
first assumption is regarding optimal doping. Rather than es¬ 
timate the doping dependence of thermoelectric performance, 
we note that in most high performance thermoelectrics, the 
Seebeck coefficient magnitude at optimal doping is between 
200 and 300 pV/K. The reasons for this are two fold. Firstly, 
Seebeck coefficients below this range do not permit high ZT. 
Note that the Wiedemann-Franz relation implies that a min¬ 
imum thermopower of 156 pV/K is necessary for a ZT of 
unity, and this is for a nil lattice thermal conductivity, which 
is clearly unrealistic. For example, the thermopower of opti¬ 
mally doped Bi 2 Te 3 is approximately 200 pV/K^. Secondly, 
for thermopower values above 300 p V/K, the chemical poten¬ 
tial is typically in the band gap, and electrical conductivity is 
correspondingly reduced due to the low carrier concentration. 
Since the electrical conductivity we measured in this material 
is already somewhat low for a high performance thermoelec¬ 
tric (about l/6th the value for optimized Bi 2 Te 3 ^), we will 
take the higher carrier concentration thermopower of 200 p 
V/ K as representing a sample of likely optimal doping. 

Our second assumption concerns the thermal conductivity. 
This may generally be written as a sum of lattice and elec¬ 
tronic thermal conductivity. The electronic thermal conduc¬ 
tivity ^electronic Can be readily estimated from the electrical 



FIG. 8. The calculated density-of-states of AgBiSe 2 . 


conductivity c and the Wiedemann-Franz relation, in which 
Keiectronic = LqOT, with Lq the Lorenz number = 2.45 x 10^^ 
(V/K)^. This relation is usually accurate for metals and heav¬ 
ily doped semiconductors as considered here. The lattice ther¬ 
mal conductivity Ktantce typically follows a 1/T relation, lead¬ 
ing to an estimation of Kiattice for AgBiSe 2 at 300 K of 0.7 
W/m-K. 

We also make a final assumption regarding the electrical 
conductivity, in particular the scattering time. The first prin¬ 
ciples calculations yield for the n-type sample experimentally 
measured a c/x at 300 K of 2 x 10^^ (f2 - m - s)^\ and given 
the measured conductivity of 160 S/cm, yield an average scat¬ 
tering time T of 8 X lO^^^s. To translate this time to a scat¬ 
tering time for optimally doped p-type requires consideration 
of two specific issues: the change in carrier type from n-type 
to p-type, and the change in carrier concentration from the 
n = 2.3 X lO^^cm^^ in our experimental work to the likely 
optimal p-type doping p = 1.2 x 10^** cm^^. 

Regarding the carrier change from n-type to p-type, one 
is really asking about the associated difference in electron- 
phonon coupling, since such coupling is the basic mechanism 
of resistivity in heavily doped semiconductors above the De¬ 
bye temperature. To help assess this issue, in Figure 8 we 
present the calculated density-of-states of AgBiSe 2 . In both 
the valence and conduction band the Se atoms contribute sub¬ 
stantially to the DOS, with the Ag comparatively less in both. 
The Bi DOS is more substantial in the conduction band than 
in the valence band. Overall, however, the relative propor¬ 
tions of the partial DOS in the conduction and valence bands 
are similar - both are substantially hybridized This suggests 
that the interatomic interactions which produce the electron- 
phonon coupling may not be too dissimilar from the conduc¬ 
tion and valence bands, and that, as a first approximation, one 
may take the electron phonon coupling and associated scat¬ 
tering times to be the same for both bands. Hence we may 
estimate the conductivity and ZT of optimally doped p-type 
AgBiSe 2 from this scattering time. We find a 300 K conduc¬ 
tivity of 800 (H-cm)^ *, approximately 80 percent of the value 
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of optimally doped 612 X 03 , which yields a 300 K ZT value of 
0.7. 

At a further level of refinement, we consider that the likely 
optimal doping level of p-type AgBiSe 2 is substantially higher 
- more than a factor of five larger - than the experimental 
carrier concentration in our n-type sample, which may affect 
carrier mobility. To assess this effect, we examine previous 
worb^ which finds that in the heavily-doped degenerate limit 
the carrier mobility p typically follows the approximate rela¬ 
tionship p ^ Taking this same proportionality for the 

scattering time yields a 300 K conductivity of 300 (fi-cm)^ * 
and an estimated ZT value of 0.4. 

We note that experimental values of ZT will depend on 
actual scattering times, so that the above should be taken only 
as first estimates. However, they demonstrate the potential 
of p-type AgBiSe 2 as a room temperature thermoelectric. 
In fact, the actual performance of this material may well be 
better than these estimates, because no effort has been made 
to model reductions of lattice thermal conductivity as are 
often attained by alloying or nanostructuring. For example, 
we may estimate the lattice thermal conductivity of optimally 
doped Bi 2 Te 3 from the Wiedemann-Franz relation, its ZT 
of unity, and Seebeck coefficient and conductivity of 200 
pF jK and^ 1000 (H-cm)^' as just 0.45 W/m-K, which is 
likely near the minimum possible for this system, and much 
lower than the 1.7 W/m-K bulk valued. We anticipate similar 
reductions may be available for AgBiSe 2 ; indeed, recent work 
on AgBiSe 2 nanoplates^ finds lattice thermal conductivity 
values of 0.45 W/m-K in the pure material and values as low 
as half this when Sb alloying is considered. 

Appendix B: Demonstration of isotropy of thermopower 
within effective mass approximation 

We here demonstrate directly that an anisotropic effective 
mass tensor leads to an isotropic Seebeck coefficient, within 
a single band model constant scattering time approximation 
(CSTA). Earlier references to this Seebeck isotropy may be 
found in Refs, and [sH We begin with the anisotropic 
dispersion relation for a single band extremum within the 
effective mass approximation (we take Ti — \ throughout); 


Ek = k^/2mx-\-k^/2my-\-k\/2m^. Now, equation (4) in Sec¬ 
tion IV yields a transport function <3aa- 

Oaa = X / '^dS (20) 

J kkl 

We now make a scale transformation, writing ka = v^2m^/4- 
and convert k'^^ky and k' to spherical coordinates. Now, the 
surface area element dS can easily be shown to be 

dS^2Efi{mx,my,m^,Q,(^) (21) 

defining an angular function fi. Similarly, one finds that 

|vk| = \/Egi{mx,my,m^,Q,^) 

defining a function gi, and 

vl,k = 2Enl{d,(^)/ma ( 22 ) 

where nc((9,‘l)) is the direction cosine of (i.e., n = cos( 0 ) 
for a = z). Now we have 

Oaa(E)= J = g2,aimx,my,m^,Q,(^)dQd(^ 

where the function g 2 ,a contains all the angular k'-space de¬ 
pendence and effective mass dependence. The angular inte¬ 
grations yield a third function gj ai'^xjiny^m^) so that we have 

Oaa{E) =E^^^g2^a{mx,my,my) (23) 

However - and this is the key point - when inserted into Eq. 
( 1 ), the effective-mass dependent term g 2 ,^a cancels, since 
it is present both in the numerator and denominator of Eq. 
(1). Hence the Seebeck coefficient is isotropic, even with an 
anisotropic effective mass tensor. 

Note that although we have assumed the constant scattering 
time approximation for simplicity, this result holds so long 
as the scattering time is a function of energy alone, and not 
direction. 
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